Magnetoresistance oscillations in multilayer systems - triple quantum wells 
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Magnetoresistance of two-dimensional electron systems with several occupied subbands oscillates 
owing to periodic modulation of the probability of intersubband transitions by the quantizing mag- 
netic field. In addition to previous investigations of these magneto-intersubband (MIS) oscillations 
in two-subband systems, we report on both experimental and theoretical studies of such a phe- 
nomenon in three-subband systems realized in triple quantum wells. We show that the presence 
of more than two subbands leads to a qualitatively different MIS oscillation picture, described as 
a superposition of several oscillating contributions. Under a continuous microwave irradiation, the 
magnetoresistance of triple- well systems exhibits an interference of MIS oscillations and microwave- 
induced resistance oscillations. The theory explaining these phenomena is presented in the general 
form, valid for an arbitrary number of subbands. A comparison of theory and experiment allows us 
to extract temperature dependence of quantum lifetime of electrons and to confirm the applicability 
of the inelastic mechanism of microwave photoresistance for the description of magnetotransport in 
multilayer systems. 



PACS numbers: 73.40.-c, 73.43.-f, 73.21.-b 

I. INTRODUCTION 

Studies of magnetoresistance, in particular, the inves- 
tigation of Shubnikov-de Haas (SdH) oscillations in semi- 
conductors and metals is an important tool for gath- 
ering information about band structure, quantum life- 
times of electrons and interaction mechanisms In 
two-dimensional (2D) electron systems, the SdH oscil- 
lations occur because of a periodic modulation of elec- 
tron scattering as the Landau levels consecutively pass 
through the Fermi level. With increasing temperature, 
when the thermal broadening of the Fermi distribution 
exceeds the cyclotron energy hu) c , the SdH oscillations are 
strongly damped. In quantum wells with at least two oc- 
cupied 2D subbands, the magnetoresistance exhibits an- 
other kind of oscillating behavior, the so-called magneto- 
intersubband (MIS) oscillations |2J. These oscillations 
occur because of a periodic modulation of the probability 
of transitions between the Landau levels belonging to dif- 
ferent subbands. The MIS oscillation peaks correspond 
to the subband alignment condition A = nfku c , where A 
is the subband separation, because the isoenergetic (elas- 
tic) scattering of electrons between the Landau levels is 
maximal under this condition. Since the origin of the 
MIS oscillations is not related to the position of Landau 
levels with respect to the Fermi energy, these oscillations 
survive at high temperatures when the SdH oscillations 
are completely damped. Early experimental studies of 
MIS oscillations have been carried out in single quan- 
tum wells with two populated 2D subbands [3]- 5]. Re- 
cently, MIS oscillations with large amplitudes have been 
observed and investigated in two-subband systems based 



on double quantum wells (DQWs) [|,0]. The DQWs ap- 
pear to be the most convenient systems for experimen- 
tal studies of this penomenon, because the two-subband 
occupation is attainable at relatively small electron den- 
sities, the electron mobility is high, and a strong tunnel 
coupling between the wells enables a high probability of 
intersubband scattering. The studies of MIS oscillations 
in DQWs provide information about temperature depen- 
dence of the quantum lifetime of electrons in the region 
where SdH oscillations are absent |6j. 

The MIS oscillations are also interesting owing to 
their interplay with another magneto-oscillatory phe- 
nomenon recently discovered in high-mobility 2D layers. 
If a 2D electron system is exposed to a continuous mi- 
crowave irradiation, microwave-induced resistance oscil- 
lations (MIROs) occur, which are governed by the ratio 
of the radiation frequency uj to the cyclotron frequency 
<jj c Q. With increasing radiation intensity, the minima 
of these oscillations evolve into "zero-resistance states" 
[3, [13] in samples with ultrahigh electron mobility. Sim- 
ilar to MIS oscillations, MIROs originate from a peri- 
odic modulation of the probability of electron transitions 
between different Landau levels. In single-subband sys- 
tems, such transitions occur because of electron scatter- 
ing in the presence of microwave excitation, when elec- 
trons absorb radiation quanta and gain the energy nec- 
essary for the transitions. A detailed theoretical descrip- 
tion of MIROs involves consideration of several micro- 
scopic mechanisms of photoresistance, which satisfactory 
describe the observed periodicity and phase of these os- 
cillations [H|-[l3|- Among them, the inelastic mecha- 
nism [13l| . associated with a microwave-generated non- 
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equilibrium oscillatory component of the isotropic part 
of the electron distribution function, dominates at low 
temperatures T, because its contribution is proportional 
to the inelastic relaxation time Ti n oc T~ 2 . Recently, 
MIROs have been studied in systems with two occupied 
subbands (DQWs) [l5j . It is found that the interplay 
between MIS oscillations and MIROs manifests itself as 
an interference of these kinds of oscillations, which is for- 
mally expressed as a product of the corresponding oscil- 
lating factors [l5[ • The observed magnetoresistance pat- 
tern strongly depends on frequency lu and exhibits inver- 
sion or enhancement of certain groups of MIS peaks. It 
is established that the inelastic mechanism of microwave 
photoresistance explains magnetoresistance oscillations 
in such two-subband systems [Hj]. 

Previous studies of MIS oscillations and their interfer- 
ence with MIROs have been restricted to two-subband 
systems. In this paper we present experimental and 
theoretical studies of these phenomena in systems with 
three occupied subbands formed in triple quantum wells 
(TQWs). The symmetric triple- well structure under in- 
vestigation is shown in Fig. [T] The barriers dividing 
the wells are thin enough to have a strong tunnel hy- 
bridization of electron states in different wells. As a re- 
sult, there exist three subbands with different quantiza- 
tion energies Ej (j = 1, 2, 3) and all of them are occupied 
by electrons at the chosen (high enough) electron density. 
We have found that the magnetoresistance of such sys- 
tems exhibits MIS oscillations with several periods deter- 
mined by subband separation energies A^y = \£j — Ej>\. 
The peculiar MIS oscillation picture is distinct from that 
observed in DQWs, where only one MIS period exists. 
This feature has not been mentioned in previous studies 
of TQWs by other groups [lj| [l7|, which concentrated 
on the regime of high magnetic fields and quantum Hall 
effect. Next, we have demonstrated that the MIS oscilla- 
tion picture in TQWs exposed to microwave irradiation 
changes and depends on the radiation frequency. This 
behavior is basically similar to that in the case of DQWs 
described above. To explain the observed MIS oscilla- 
tions and their interference with MIROs in TQWs, we 
generalize the magnetoresistance theory in the presence 
of microwave irradiation to the multisubband case. All 
experimental results are in agreement with theoretical 
caluclations involving the inelastic mechanism of pho- 
toresistance. 

The paper is organized as follows. Section II presents 
our experimental results for TQWs and a theoretical de- 
scription of the MIS oscillations in many-subband sys- 
tems. Section III presents the magnetoresistance under 
microwave irradiation, also together with a theoretical 
description, a comparison of theory and experiment, and 
a discussion of the results. Conclusions are given in the 
last section. Appendix A contains details of the theo- 
retical calculation of photoresistance for many-subband 
systems, and Appendix B describes the tight-binding ap- 
proach for a calculation of subband spectrum and scat- 
tering rates in symmetric triple- well systems, with appli- 
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FIG. 1: (color online) (a) Image of a triple quantum well and 
(b) Landau level staircase for a triple quantum well with three 
occupied subbands (1,2,3). 



cation to our samples. 



II. MIS OSCILLATIONS IN TRIPLE QUANTUM 
WELLS 

Our samples are symmetrically doped GaAs TQWs, 
separated by Al^Gai-^As barriers, with a high total elec- 
tron sheet density of n s — 9 x 10 11 cm~ 2 and mobilities 
of 5 x 10 5 cm 2 /V s (wafer A) and 4 x 10 5 cm 2 /V s (wafer 
B). The central well width is about 230 A and both side 
wells have equal widths of 100 A. The barrier thickness 
db is 14 A (wafer A) and 20 A (wafer B). In order to make 
the central well populated, we increased its width. The 
estimated density in the central well is 35% smaller than 
in the side wells. The layers are shunted by ohmic con- 
tacts. Figure [1] gives an image of a triple quantum well 
with three occupied subbands (j = 1, 2, 3) and the stair- 
cases of Landau levels. The subband separation energies 
Ajji , which characterize the coupling strength between 
the quantum wells (see Appendix B), are presented in 
Table HI The measurements have been carried out in a 
dilution refrigerator (low temperature MIS studies) and 
in a VTI cryostat using a waveguide for microwave ex- 
periments to deliver microwave radiation down to the 
sample. A conventional lock-in technique for magneto- 
transport measurements under a continuous microwave 
irradiation (35 GHz to 170 GHz) has been used. Several 
specimens of both van der Pauw and Hall bar geometries 
from both wafers have been studied. 



TABLE I: Subband separation energies for two samples, ex- 
tracted from Fourier analysis of magnetoresistance. 
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The description of the MIS oscillations will be focused 
on samples with dh — 14 A owing to a stronger tunnel 
coupling which gives rise to better pronounced MIS fea- 
tures. In Fig. [5] (a) we present temperature dependence 
of MIS oscillations from T = 1.4 K to T = 4.2 K. The 
inset also shows MIS oscillations at T = 50 mK. For 
this low temperature and also at 1.4 K, the MIS oscil- 
lations are superimposed on low-held SdH oscillations. 
Figure [2] (b) demonstrates MIS oscillations for both sam- 
ples with df, = 14 A and df, = 20 A . It is obvious that for 
the thicker barrier, when tunnel coupling is weaker, the 
subband separation A^y becomes smaller and the prob- 
ability of intersubband transitions of electrons decreases. 
This is reflected in the periodicity and in the amplitude 
of MIS oscillations (see also Table A} . 
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FIG. 2: (color online) (a) Temperature dependence for MIS 
oscillations in a TQW with d b = 14 A for 1.4, 2.0, 3.2 and 
4.2 K. The inset shows MIS oscillations at 50 mK, super- 
imposed on low-field SdH oscillations, (b) Comparison of 
both wafers with db = 20 A (top, shifted up for clarity) and 
d b = 14 A (bottom) at T = 4.2 K. 

To describe the data in more detail, we have general- 
ized the theory of magnetoresistance in the systems with 
two occupied subbands (DQWs) jf| [H, [li| to the case of 
JV subbands. We consider elastic scattering of electrons 
in the presence of a magnetic field under the condition 
of large filling factors (Fermi energy ep is much larger 
than hui c ), and apply the self-consistent Born approxi- 
mation to describe the density of states and the linear 
response. The expression for magnetoresistance is conve- 
niently presented in the form 



(o) 
Pd = Pd 



(2) 
Pd > 
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where p d °^ is the classical resistivity, pp is the first-order 
(linear in Dingle factors) quantum contribution describ- 

( 2) 

ing the SdH oscillations, and p d ' is the second-order 



(quadratic in Dingle factors) quantum contribution con- 
taining the MIS oscillations. In the regime of classically 
strong magnetic fields, one obtains 
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where m is the effective mass of electrons, dj = 
cxp(— nvj/Lo c ) are the Dingle factors, T = X/sinhX 
with X — 2TT 2 T/huj c is the thermal suppression fac- 
tor, and rij are the partial densities in the subbands 
GZj n j ~ n s)- The sums are taken over all subbands, 
and since for the terms with j = j' one has A.jji = 0, 
the corresponding cosines in Eq. (4) are equal to 1 . The 
subband-dependent quantum relaxation rates Vj and v^y , 
as well as the transport scattering rates v 1 -" and v% en- 
tering Eqs. (2)-(4) are defined according to 
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where Wjj> (q) are the Fourier transforms of the cor- 

1 - 

2kjkji cos 0/(kj + k 2 , ), 9 is the scattering angle, and 



relators of the scattering potential, Fjji(0) 
9/(k 2 + k 2 ,), 9 is the scattei 

k 3 = ^27 



inn. 



is the Fermi wavenumber for subband j. 
Since Fjj = 1 — cos 9, the intrasubband transport rates 
Ujj are defined in a conventional way. The theory is valid 
if the subband separations Ajji are large compared to the 
broadening energies hvj and the Dingle factors are small, 
d% <C 1. In a similar way as we introduced the averaged 
transport time r tr by Eq. (2), one can introduce the av- 
eraged quantum lifetime l/r q = N~ x Y] ^ Vj. Application 
of Eqs. (l)-(6) to the particular case of three subbands 
(JV = 3, j = 1, 2, 3) is straightforward. 

The behavior of magnetoresistance can be illustrated 
within a simple model using equal electron densities rij = 
n s /JV and assuming that all and dj are equal to each 
other, in particular, dj — d = exp(—Tr/uj c T q ). Neglecting 
the SdH oscillations, we obtain for N = 3 
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The MIS oscillations are represented as a superposition 
of three oscillating terms determined by relative positions 
of the subband energies. Notice that this expression does 
not depend on transport rates except the one standing 
in the Dingle factor d. The approximation (7), in prin- 
ciple, can be applied for estimates to our system, since 
we have high total electron-sheet density and a strong 
tunnel coupling. To describe experimental magnetoresis- 
tance in detail, a more careful calculation based on Eqs. 
(l)-(6) is required. 

We calculate the magntoresistance of our system under 
a simplified assumption that the scattering potential is 
essential only in the side (s) wells, since the growth tech- 
nology implies that most of the scatterers reside in the 
outer barriers. The correlation length l c = 18.3 nm en- 
tering the scattering potential correlator Wjj>(q) ex w s (q) 
(see Appendix B for details) is determined by comparing 
the results of calculations to low-temperature magnetore- 
sistance data for the samples with mobility 5 x 10 5 cm 2 /V 
s. The averaged quantum lifetime estimated in this way 
is T q ~ 3.8 ps. The experiment shows a slow suppression 
of the MIS oscillations with temperature, which occurs 
owing to the contribution of electron-electron scattering 
into Landau level broadening. Though the theory pre- 
sented above does not take this effect into account ex- 
plicitly, it can be improved by replacing the quantum 
relaxation rates according to 

T 2 

Vj -> Vj + v ee , v ee = A- — (8) 

hep 

where v ee is the electron-electron scattering rate 0, [2l[ , 
the Fermi energy is expressed through the averaged elec- 
tron density as Ef = h 2 ir(n s /3)/m, and A is a numerical 
constant of order unity. In Fig. [3] we present a compar- 
ison of experiment and theory for two choosen tempera- 
tures, T = 6 K and T = 10 K. We have carried out this 
procedure for many temperatures from T = 1 K up to 
30 K, and estimated v ee by fitting the amplitudes of the- 
oretical and experimental magnetoresistance traces. The 
effect of electron-electron scattering becomes essential for 
T > 2 K and strongly reduces the amplitude of the MIS 
oscillations at T ~ 10 K. As seen from the log-log plot 
in the inset to Fig. [3l the extracted scattering rate v ee 
follows the T 2 dependence, in accordance with Eq. I 
This behavior is similar to that observed in DQWs 
Using ef — 10.5 meV, we find A = 2.2. 



III. INFLUENCE OF MICROWAVES ON MIS 
OSCILLATIONS 

In this section, we investigate the influence of a con- 
tinuous microwave irradiation on our TQW system. We 
have studied power, temperature, and frequency depen- 
dence of magnetoresistance for both wafers, though we 
focus again on the samples with = 14 A . 

In the upper part of Fig. [T] we present the magnetore- 
sistance for different microwave powers at a temperature 
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FIG. 3: (color online) Comparison of the experimental and 
theoretical traces for a TQW with d b = 14 A at T = 6 K 
(top) andT = 10 K (bottom). By fitting the amplitude of MIS 
oscillations, electron-electron scattering rate u e( , is extracted 
(see the points in the inset). The linear fit to experimental 
data (line) corresponds to the theoretical dependence of Eq. 
(8) with A = 2.2. 



of T = 4.2 K and a frequency of 55 GHz. Without mi- 
crowave irradiation (no MW), only the MIS oscillations 
are visible. An increase in microwave power (—10 dB 
attenuation) leads to an enhancement of all MIS fea- 
tures for B < 0.25 T and to a damping of all such 
features for B > 0.35 T whereas the MIS oscillations 
around B = 0.3 T are almost unchanged. A further 
increase in power (—2 dB attenuation) leads to a damp- 
ing of the MIS oscillation amplitude for B < 0.25 T, 
slightly increased compared to the MIS oscillation ampli- 
tude without microwave irradiation. The MIS features 
around B = 0.3 T are considerably damped, while for 
0.35 T < B < 0.5 T the MIS peaks are inverted. No 
polarization dependence of magnetoresistance has been 
found. 

A similar behavior, with enhanced, suppressed, or in- 
verted MIS peaks is observed in the magnetoresistance 
measured at different microwave frequencies in the range 
between 35 GHz and 170 GHz, as shown in the up- 
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FIG. 4: (color online) Measured (upper panel) and calculated 
(lower panel) magnetoresistance for TQW with db = 14 A at 
a lattice temperature of T = 4.2 K under excitation by mi- 
crowaves with a frequency of 55 GHz. The microwave electric 
field used for calculation of the magnetoresistance (see de- 
tails in the text) corresponds to a different microwave power 
in dB. 



per part of Fig. [5] A strongly modified picture of the 
Mf S oscillations correlates with the microwave frequency. 
The features most affected by the microwave irradiation, 
strongly sensitive to its frequency occur at B = 0.27 T 
and B = 0.43 T. The plots for 110 GHz and 170 GHz 
definitely show several regions of enhanced peaks and 
two regions of suppressed or inverted peaks (for exam- 
ple, the regions around 0.18 T and 0.34 T for 170 GHz). 
For 35 GHz, all the MIS oscillations above 0.2 T are in- 
verted. For 35 GHz and 70 GHz, some SdH oscillations 
are visible in the region above 0.4 T. For 110 GHz and 
170 GHz, when the absorption of microwave radiation in 
this region is higher, the SdH oscillations are suppressed 
because of the heating of the electron gas by microwaves. 

The peculiar features of the microwave-modified mag- 
netoresistance can be understood in terms of interference 
of the MIS oscillations with MIROs. This phenomenon is 
already known for two-subband systems [15j . The theory 
applied for explanation of our measurements is presented 
below. The dissipative resistivity in the presence of mi- 



FIG. 5: (color online) Measured (upper panel) and calculated 
(lower panel) magnetoresistance for a TQW with db = 14 A 
at the lattice temperature T — 1.4 K for several choosen fre- 
quencies. The curves for 70, 110, and 170 GHz are shifted up 
for clarity. 



crowaves is given by (see Appendix A for derivation of 
the expressions) 



Pd = Pd + Pin + Pdi, 



(9) 



where the dark resistivity pd is described in Sec. II, while 
Pi n and pdi are the microwave-induced contributions due 
to inelastic and displacement mechanisms, respectively. 
At low temperatures, pi n is the main contribution. It is 
given by 

vjvjdjdji cos - — , (10) 
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is an oscillating function describing MIROs, and 
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is the dimensionless function proportional to the mi- 
crowave power. In this expression, v F = N~ z2j Vj is 
the averaged Fermi velocity, and Ti n is the averaged in- 
elastic relaxation time introduced in Appendix A. The 
averaged transport time r tr is defined in Sec. II, and the 
coefficients s± are given in Appendix A. Notice that in 
the case of linear polarization of radiation, and away from 
the cyclotron resonance, |s+| 2 + |s~| 2 — {uj 2 + ui 2 ) / (uj 2 — 
uj 2 ) 2 . The contribution p^ is presented in Appendix A. 
Since pdi is much smaller than pi n , it is not taken into 
account in our consideration. The influence of microwave 
radiation on the SdH oscillations can also be neglected 
at weak radiation power. 

To get a visual description of the influence of radia- 
tion on magnetoresistance, we again use a simple model 
assuming equal partial densities rij — n s /N , equal trans- 
port scattering rates u% and Dingle factors dj = d. The 
magnetoresistance of three-subband system (N = 3) then 
takes the form 
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which differs from Eq. (7) only by the presence of the fac- 
tor 1—A U . The products of A u by the MIS oscillation fac- 
tors cos(2tv Ajf /Hue) lead to interference oscillations of 
the magnetoresistance. In the regions of frequency where 
A u is negative, one expects an enhancement of the MIS 
peaks. If A^ is positive, the peaks are suppressed and 
inverted with increasing microwave power. This is the 
main feature of the behavior we observe experimentally 
in Figs. [landEl 

A comparison of experimental results with theory 
based on Eqs. (10)-(12) is demonstrated in the lower 
parts of Figs. [4] and [5] Apart from the known parame- 
ters used also in Sec. II, we apply the following estimate 
for the inelastic relaxation time [l3j]: t„ ~ Hef/T 2 , as- 
suming that the relaxation is governed by the electron- 
electron interaction. The reliability of this estimate 
is confirmed in numerous experiments on magnetoresis- 
tance influenced by either microwave field [H, [22[ or 
static electric field [H, 24]. To explain the experimen- 
tal data, it is important to take into account microwave 
heating of the electron gas. This effect is directly visible 
in our experiment and results in a suppression of the SdH 
oscillations under microwave irradiation. The increase of 
the effective electron temperature over the lattice tem- 
perature also leads to a decrease of the inelastic relax- 
ation time and quantum lifetimes, see Eq. (8), so the 
MIS oscillation amplitudes are expected to be suppressed 
as a result of electron heating. The electron temperature, 
which depends on the magnetic field, radiation frequency, 
and power, has been calculated assuming energy relax- 
ation of electrons due to their interaction with acoustic 
phonons. Finally, to determine the electric field E u cor- 
responding to our measurements, we use an estimate for 



the microwave electric field generated by our source as 
10 V/cm (at 55 GHz). Thus, the attenuations of -10 dB 
and —2 dB correspond to E u — 3.2 V/cm and 7.9 V/cm, 
respectively, and we applied these values for calculation 
of the magnetoresistance shown in Fig. [4] 

The theoretical plots in Fig. 0] reproduce all the basic 
features of the experimental magnetoresistance traces, in 
particular, a suppression and inversion of two MIS peaks 
around B = 0.4 T, because of the contribution pi n with 
positive A w . Notice that the non-monotonic power de- 
pendence of the MIS peaks around 0.2 T is explained by 
the interplay of MIS /MIRO interference and heating ef- 
fects. At low radiation power the enhancement of these 
peaks occurs because of the contribution p in with nega- 
tive A,,,. At high power, when the saturation effect takes 

place mm* 

a decrease in the Dingle factors due to the 
heating-induced increase in v ee becomes more important 
and the MIS peaks are suppressed. 

The expected microwave electric field in the frequency- 
dependent measurements shown in Fig. [5] is E u ~ 3 
V/cm. To get a closer resemblance of the theoretical 
magnetoresistance to the experimental plots, we slightly 
varied E u around this value and obtained the best fit at 
E u = 3.5 V/cm for 35 GHz and 70 GHz, 4 V/cm for 
110 GHz, and 2.2 V/cm for 170 GHz. The correspond- 
ing theoretical plots are presented in Fig. \5\ Since the 
lattice temperature for these measurements is 1.4 K, the 
heating effect appears to be considerable. For 35, 70, 
and 110 GHz, the calculated electron temperature in the 
vicinity of the cyclotron resonance is about 3.5 K, which 
is close to our experimental estimates obtained from sup- 
pression of the SdH oscillations. In general, a reasonably 
good agreement between theory and experiment at dif- 
ferent frequencies suggests that the theoretical model ap- 
plied for the calculations is reliable. 



IV. CONCLUSIONS 

We have studied transport properties, including mi- 
crowave photoresistance, of the electron systems with 
three occupied 2D subbands in perpendicular magnetic 
fields. Such systems are realized in TQWs with high 
enough electron density. As we have demonstrated, both 
experimentally and theoretically, the magnetoresistance 
of TQWs is qualitatively different from that for single- 
subband and two-subband systems, and contains a super- 
position of three oscillating terms whose frequencies are 
given by the subband separation energies A12, A 13, and 
A23. This occurs because the quantum contribution to 
the resistivity is essentially determined by electron scat- 
tering between the Landau levels of different subbands. 
Therefore, there exist MIS oscillations of resistivity, and 
the picture of these oscillations becomes complicated in 
the systems with more than two occupied subbands. We 
have presented a theoretical description of such oscilla- 
tions by generalizing the theory of quantum magnetore- 
sistance to the multisubband case, and obtained a good 
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agreement with the experiment. 

Similar as in single-subband and two-subband sys- 
tems, the quantum contribution to the resistivity de- 
creases with increasing temperature T because of the 
decrease in quantum lifetime due to enhanced contribu- 
tion of electron-electron scattering. By measuring the 
amplitude of the MIS oscillations at different temper- 
atures up to 30 K, we have established that the tem- 
perature dependence of electron-electron scattering rate 
follows the theoretically predicted T 2 law, see Eq. (8). 
The numerical constant A in this dependence, A = 2.2, 
is close enough to those determined from the MIS os- 
cillations in two-subband systems, both in double quan- 
tum wells (A = 3.5) and in single quantum wells (2f| 
(A = 2.6). Therefore, one may conclude that the influ- 
ence of electron-electron scattering on the quantum life- 
time of electrons is not very sensitive to the number of 
occupied subbands. 

The TQWs exposed to a continuous microwave irradia- 
tion demonstrate dramatic changes in magnetoresistance. 
The effect of microwaves is understood as a result of in- 
terference of MIS oscillations and microwave-induced re- 
sistance oscillations. A similar effect takes place for two- 
subband systems in double quantum wells 15j , where it is 
easier recognizable owing to a simpler picture of the MIS 
oscillations. To describe our observations, we have de- 
veloped a theory of magnetoresistance of multisubband 
systems under microwave irradiation, and applied it to 
our three-subband systems. Among the mechanisms of 
microwave photoresistance, the inelastic mechanism is 
found to be responsible for the observed magnetoresis- 
tance features. In spite of several approximations of the 
theory, in particular, those for description of elastic scat- 
tering (see Appendix B) , we have obtained a good agree- 
ment between theoretical and experimental magnetore- 
sistance traces by using an established estimate for the 
inelastic relaxation time. 

In summary, our investigation of low-field magneto- 
transport in three-subband systems both with and with- 
out microwave excitation is a useful step towards un- 
derstanding the influence of energy spectrum and scat- 
tering mechanisms on the transport properties of low- 
dimensional electrons. 
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APPENDIX A: MICROWAVE 
PHOTORESISTANCE OF A MANY-SUBBAND 
SYSTEM 

In the presence of electromagnetic radiation (mi- 
crowaves) of frequency to and under a dc excitation, one 
can derive the quantum Boltzmann equation for electrons 



in a magnetic field by using a transition to the moving co- 
ordinate frame, in a similar way as for the single-subband 
system (see [lj] and references therein). This leads to 
the kinetic equation for the Wigner distribution function 
fjetp, which depends on the subband index j, energy e, 
and angle ip of the electron momentum: 



dip 



E 

i 



+nLU+y jj , ,ip' fjeip] + ^m-(Al) 

Notice that in this Appendix we use the system of units 
where h = 1. In the kinetic equation, we introduced the 
dimensionless (normalized to its zero field value) density 
of states Dj(e). Next, J n (x) is the Bessel function, J in is 
the collision integral describing inelastic scattering, and 
Vjji are the scattering rates defined in Sec. II. The other 
quantities standing in Eq. (Al) are 



+ s+(vje~ l<p 



Vj>e 



(A2) 



and 



■E+fae-^-Vfe-**)], (A3) 



where E± 



E x ± iE y , E 



(E x ,E y ) is the dc field 



strength, is the strength of microwave electric field 
(related to the incident microwave field strength Ei in 
vacuum as E u = Ei/\/e*, see below), and Vj — pj/m 
are the subband-dependent Fermi velocities. The factors 
s± describe polarization of the radiation and account for 
electrodynamic effects (2(| [2?]]. For the case of linear 
polarization, 



s± 



(A4) 



where lu p — 2ire 2 n s /mc\i / €* is the plasma frequency, 

= (V e «ac + V<^)/ 2 ' £ vac = 1 is the dielectric per- 
mittivity of vacuum, and is the dielectric permittivity 
of the medium surrounding the quantum wells. 

The distribution function can be expanded in the an- 
gular harmonics: fj etp = J^i fjsis ll(p ■ The density of dis- 
sipative electric current in the 2D plane, j = (j x ,jy), is 
determined by the I = 1 harmonic: 



: 3x 1]y 



j * deDj(e)f jel , (A5) 



In the regime of classically strong magnetic fields, the 
anisotropic part of the distribution function is expressed 
through the isotropic (angular-independent) part fj £ = 
fj e i for I — 0. This leads to the expression for the current 
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in the form 



jj' 



i ; 27T Jo 2 7T 



-/is].(A6) 

The response to E— \j- = OdE-] gives the symmetric 
part of dissipative conductivity considered below. The 
resistivity is then given by p h d IW = (Td/cr\, where a± = 
e 2 n s /mu! c is the classical Hall conductivity. 

The isotropic part of the distribution function can be 
represented in the form 



fie — f t 



(0) 



(A7) 



where fe is a slowly varying function of energy, which 
is close to a quasi-equilibrium (heated Fermi) distri- 
bution, while gj E is a rapidly oscillating (periodic in 
Hl> c ) function, which is also represented as gj £ — 
Sfe 9jk exp(27rifce/w c ). After a substitution of expres- 
sion (A7) into Eq. (A6), the contribution coming from 

fe ' 1 produces the "dark" resistivity and its modification 
by the microwaves due to displacement mechanism. The 
microwave modification of the resistivity due to the in- 
elastic mechanism originates from the term proportional 
to gj e . In the following, we search for the response linear 
in the dc field and quadratic in the microwave field. This 
is done by expanding the Bessel functions in powers of 
(3jf and retaining only the lowest-order terms. Also, we 
use the lowest-order expansion of the density of states in 
the Dingle factors, Dj(e) ~ 1 — 2dj cos[27r(e — ej)/ui c ], 
so only the lowest oscillatory harmonics (k = ±1) are 
relevant. The angular and energy averaging in Eq. (A6) 
in this case are carried out analytically, and one gets the 
expression for the dark resistivity pd (Sec. II), as well as 
the microwave- induced contributions pdi and pi n : 



Pdi 



m 



e 2 n s 



sin 1 sin 

U) c uj c uj c 

° 2nA jf 

Vjjrdjdjr cos — . 



(A8) 



where P w is defined by Eq. (12), 



dO 
o 2^ 



u jr (e)[F 3 f(0)} 2 , 



and 



m 2iru> sr-^ n.j + ny tr 



E 



e n s ui c * — ' n s 

33 

( 2ui£j , 

.'/;:''; exp I— — I <:(: 



(A10) 



One can find gj e from the isotropic part of Eq. (Al) by 
using the relaxation time approximation for the isotropic 
part of the inelastic collision integral: 



J;, 



fje ft 



(0) 



(All) 



p(0) 



This approximation is valid for small deviations fj E ~fe 
and is justified in a similar way as for the single-subband 
systems, based on a linearization of the collision integral 
for electron-electron scattering [l3[ . After substitution of 
expression (All) into Eq. (Al), one can obtain a system 
of linear equations for gjk, which is easily solved under 
a reasonable condition that the intersubband scattering 
dominates over the inelastic one: Vjji ^> 1/rj" (j =^ j'). 
This gives subband-independent harmonics gjk — gk, in 
particular, 



9i = 



P u r m sin(27rw/o; c ) 
1 + P u (n n /T t r) sin 2 (7rw/w c ) 



where the averaged inelastic relaxation time is defined 
according to l/r in — N^ 1 ^^ 1/rj™. The denominator 
in Eq. (A12) describes the effect of saturation. A sub- 
stitution of the result Eq. (A12) into Eq. (A10) leads to 
Eq. (10). 



APPENDIX B: ELECTRON SPECTRUM AND 
SCATTERING RATES IN A TRIPLE- WELL 
SYSTEM 

To describe the scattering rates, we employ the wave 
functions and electron energies in the subbands found 
from the tight-binding Hamiltonian [28[ using the expan- 
sion of the wave function ip{ z ) = J^i fiFiiz) in the basis 
of single-well orbitals Fi(z) (i — 1, 2, 3 numbers the left, 
central, and right well, respectively). This leads to the 
matrix equation for the coefficients iff. 



i Ej — e —t\2 

(0) 



V 



-*12 




eh 



"^23 




(Bl) 



where s\ are the single-well quantization energies and 
(A9) ta' are the tunneling amplitudes. For the case of sym 



(0) 



(0) _ 



-(0) _ 



metric TQWs [e\ 
the energies of the three subbands, e 



£cj ti2 — *23 = t) 

are SI 



si = 
£3 = 



•e.)/2-A, 
e s )/2 + A, 



(B2) 



E s ) 2 /A + 2t 2 . 
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The corresponding eigenstates are expressed through the 
single- well orbitals as tpj(z) = J2iXijFi( z )- The matrix 
Xij is given by 

/ Cit/(e- - ei) VV2 C 3 t/(e s -e 3 )\ 
Xij = Ci C 3 , (B3) 

V Cit/(e. - ei) -1/V2 C 3 t/(e s - e 3 ) J 

where Ci, 3 = [l + 2t 2 /(e s - £i, 3 ) 2 ] . This matrix is 
composed from the three columns of ipi for the states 
J = 1,2,3. 

The parameters of the tight-binding model can be ex- 
tracted from the subband gaps found experimentally. By 
setting e s as the reference energy, one has: 

e c = A 23 -A 12 , t = v /A 23 A 12 /2. (B4) 

For our samples (wafer A) we obtain s c = 2.5 meV and 
2t = 3.35 meV. Using the total density n s = 9 x 10 11 
cm~ 2 , one can find the subband densities n\ = 3.62 x 10 11 
cm" 2 , n 2 = 3.23 x 10 11 cm" 2 , and n 3 = 2.14 x 10 11 
cm -2 . The electron density in each side well is n S id e — 
J2j Xij n j — J2j xij n j = 3-23 x 10 11 cm -2 , and the elec- 
tron density in the central well is n cent — J2j x\j n j = 
2.53 x 10 11 cm" 2 . 

The random scattering potential acting on electrons is 
V(r,z), where r is the in-plane coordinate vector, and 
z is the coordinate in the growth direction. The matrix 
elements of this potential, V^y(r), are expressed through 
the effective 2D potentials in the layers, introduced as 
Vi(r) = J dz\Fi(z)\ 2 V(r, z). Accordingly, the correlators 
of the potentials are written through the correlators of 
^(r): 

W of {\v-v'\)^{{V^)V f3 (v'))) 

= ^ J XiiXi3'Xi'j'Xi'3 w n'(\ T - r 'D- ( B5 ) 

ii' 



The product of the factors \ determines the overlap of the 
electron wave functions, while the factors Wi»'(|r — r'|) = 
((Vi(r)Vi' (r'))) describe intralayer i = i' and interlayer 
i ^ i' potential correlations. In TQWs one can neglect 
interlayer correlations between the potentials of the side 
wells because of a large distance between these wells: 
Wis ~ 0. Then, symmetric TQWs are characterized by 
three correlators: W s — Wu — W 33 , W c — W22, and 
W sc = W12 = TT 23 = TT 2 i = TT 32 . A reasonable approx- 
imation for our TQWs is to assume that the effective 
potential in the central well, V 2 (r), is much weaker than 
the side-well potentials. In this case, W c and W sc can be 
neglected compared to W s , and Eq. (B5) is rewritten as 



TT jy (|r-r'|)^2 X 2 jX 2 J ,iy s (|r-r'|). (B6) 



The spatial Fourier transform Wjji (q) = 
Jdrexp (— iq • r) Wjj/(|r|) determines all scattering 
rates according to Eqs. (5) and (6). In the approxi- 
mation Eq. (B6), Wjj>(q) ~ xijXij'W s {q), where w s (q) 
is the spatial Fourier transform of W s . The concrete 
form of the function w s (q) depends on the nature of 
the scatterers, their distribution in the structure, and 
on the TQW potential which determines the shape of 
Fi(z). In the case of long-range scattering potential, 
the magnetoresistance pd(B)/p d (0) is weakly sensitive 
to this form, but essentially depends on the effective 
correlation length l c which defines the scale of the 
g-dependencc. In our calculations, we use a model 
w s (q) cx exp(-l c q). 
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